Effective low-energy theory of superconductivity in carbon nanotube ropes 
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We derive and analyze the low-energy theory of superconductivity in carbon nanotube ropes. 
A rope is modelled as an array of metallic nanotubes, taking into account phonon-mediated as 
well as Coulomb interactions, and arbitrary Cooper pair hopping amplitudes (Josephson couplings) 
between different tubes. We use a systematic cumulant expansion to construct the Ginzburg-Landau 
action including quantum fluctuations. The regime of validity is carefully established, and the effect 
of phase slips is assessed. Quantum phase slips are shown to cause a depression of the critical 
temperature T c below the mean- field value, and a temperature-dependent resistance below T c . We 
compare our theoretical results to recent experimental data of Kasumov et al. [Phys. Rev. B 68, 
214521 (2003)] for the sub-T c resistance, and find good agreement with only one free fit parameter. 
Ropes of nanotubes therefore represent superconductors in the one- dimensional few-channel limit. 

PACS numbers: 73.63.Fg, 74.78.Na, 74.25.Fy 



I. INTRODUCTION 

Over the past decade, the unique mechanical, electri- 
cal, and optical properties of carbon nanotubes, including 
the potential for useful technological applications, have 
created a lot of excitement yT, U While many of these 
properties are well understood by now, the experimental 
observation of intrinsic 0- 13 and anomalously strong 
proximity-induced 0, Q superconductivity continues to 
pose open questions to theoretical understanding. In this 
paper we present a theory of one-dimensional (ID) super- 
conductivity as found in ropes of carbon nanotubes 0, Q 
and potentially in other nanowires. Ropes are ID ma- 
terials in the sense that there is only a relatively small 
number of propagating channels (typically, N sw 10 to 
100) available to electronic transport. While most other 
ID materials tend to become insulating at low tempera- 
tures due to the Peierls transition or as a consequence of 
electron-electron interactions, nanotubes can stay metal- 
lic down to very low temperatures yj. If the repulsive 
electron-electron interactions can be overcome by attrac- 
tive phonon-mediated interactions, ropes of nanotubes 
can then exhibit a superconducting transition. 

However, due to strong ID fluctuations, this transi- 
tion is presumably rather broad, and the question of how 
precisely superconductivity breaks down as the number 
of propagating channels decreases has to be answered 
by theory. Experimentally, the breakdown of supercon- 
ductivity manifests itself as a temperature-dependent re- 
sistance below the transition temperature T c , which be- 
comes more and more pronounced as the rope gets thin- 
ner Q . According to our theory, this resistance is caused 
by quantum phase slips, and therefore the experimental 
data published in Ref. [j] have in fact explored a regime 
of ID superconductivity with clear evidence for quan- 
tum phase slip events that had not been reached before. 
To the best of our knowledge, nanotube ropes represent 
wires with the smallest number of propagating channels 
showing intrinsic superconductivity, even when compared 
to the amorphous MoGe wires of diameter w 10 nm stud- 



ied in Ref. [8|, where still several thousand channels are 
available. 

We theoretically analyze superconductivity in nan- 
otube ropes by starting from the microscopic model of 
an array of N individual metallic single-wall nanotubes 
(SWNTs) without disorder, with effectively attractive 
on-tube interactions and inter-tube Josephson couplings. 
A similar model has been suggested by Gonzalez [3, E| • 
In the absence of the Josephson couplings, each SWNT 
would then correspond to a Luttinger liquid with inter- 
action parameter g c+ > 1, where g c+ = 1 marks the non- 
interacting limit. For simplicity, we take the same g c + on 
each SWNT. For example, for (10, 10) armchair SWNTs, 
assuming good screening of the repulsive Coulomb in- 
teractions, phonon exchange via a breathing mode (as 
well as optical phonon modes) leads to g c + « 1.3, see 
Ref. yd]]. In the case of attractive interactions, the dom- 
inant coupling mechanism between different SWNTs is 
then given by Cooper pair hopping, while single-particle 
hopping is drastically suppressed by momentum conser- 
vation arguments |3, Il2| . The coupling among different 
SWNTs is thus encoded in a Josephson coupling matrix 
Aij, where i,j = 1, . . . , N. As different nanotube chiral- 
ities are randomly distributed in a rope, only 1/3 of the 
SWNTs can be expected to be metallic. In general, the 
matrix should therefore be drawn from an appropri- 
ate random distribution. We consider below one individ- 
ual rope with a fixed (but unspecified) matrix, and derive 
general statements valid for arbitrary Aij. In that sense, 
our theory allows to capture some disorder effects, at 
least qualitatively. However, since typical elastic mean 
free paths in SWNTs exceed l^m yj, disorder effects 
within individual SWNTs are ignored completely. The 
above reasoning leads us to the problem of N coupled 
strongly correlated Luttinger liquids, where the number 
of "active" chains N < 100 with reference to the experi- 
ments of Ref. Q . This is a difficult problem that neither 
permits the use of classical Ginzburg-Landau (GL) the- 
ory nor of the standard BCS approach, in contrast to 
the situation encountered in, e.g., wide quasi-lD organic 
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superconductors fl3| . 

The approach taken in this paper is sketched next. Af- 
ter a careful derivation of the coupled-chain action in 
Sec. [n] we proceed by introducing the appropriate or- 
der parameter field. In Sec. IIIII we then perform a cu- 
mulant expansion in this order parameter, and thereby 
give a microscopic derivation of the quantum GL action, 
which then allows to make further progress. We estab- 
lish the temperature regime where this theory is reliable, 
and then focus on the important phase fluctuations of 
the order parameter field. At temperatures T well be- 
low a mean-field transition temperature T c °, amplitude 
fluctuations are shown to be massive, and hence the am- 
plitude can safely be treated in mean-field theory. The 
massless phase fluctuations then capture the important 
physics, and we specify the resulting effective low-energy 
action, valid at temperatures well below T®. Based on 
this action. Sec. IIVI explains why quantum phase slips 

(qpss) [M E m G3 

are crucial for an understand- 
ing of the experimental results of Refs. 0,0]. First, they 
cause a depression of the transition temperature T c below 
the mean- field critical temperature T®. Furthermore, for 
T < T c , a finite resistance R(T) due to QPSs appears, 
which exhibits approximate power-law scaling. We de- 
termine the full temperature dependence of R(T < T c ) 
for arbitrary rope length in Sec. [V] In Sec. IVII we then 
compare these results for R(T) to the experimental data 
of Ref. 3, focussing on two of their samples. Finally, 
Sec. I Vlll offers some concluding remarks. Throughout 
the paper, we put h — ks = 1. 

II. MODEL AND ORDER PARAMETER 

We consider a rope consisting of N metallic SWNTs 
participating in superconductivity. Experimentally, this 
number can be found from the residual resistance mea- 
sured as offset in the resistance when extrapolating down 
to T = 0. Due to the attached normal electrodes in 
any two-terminal measurement of the rope, despite of the 
presence of superconductivity, there will always be a fi- 
nite contact resistance R c . Since each metallic tube con- 
tributes two conduction channels, assuming good trans- 
parency for the contacts between metallic tubes and the 
electrodes, this is given by 

Rc= w> Rq = h / 2e2 - 12 - 9kQ - (2- 1 ) 

Extrapolation of experimental data for the resistance 
R(T) down to T — > within the superconducting regime 
then allows to measure R c , and hence N. Good trans- 
parency of the contacts is warranted by the sputtering 
technique used to fabricate and contact the suspended 
rope samples in the experiments of Refs. 0,0]. An al- 
ternative way to estimate N comes from atomic force 
microscopy, which allows to measure the apparent radius 
of the rope, and hence yields an estimate for the total 
number of tubes in the rope. On average, 1/3 of the 



tubes are metallic and one should obtain the same 
number N from this approach. Fortunately, these two 
ways of estimating N provide consistent results in most 
samples 0. Therefore the values for N used below are 
expected to be reliable. 

Here we always assume that phonon exchange leads 
to attractive interactions overcoming the (screened) 
Coulomb interactions. This assumption can be prob- 
lematic in ultrathin ropes, where practically no screen- 
ing arises unless there are close-by gate electrodes. For 
sufficiently large rope radius, however, theoretical argu- 
ments supporting this scenario have been provided in 
Ref. I n the absence of intra-tube disorder, then the 
appropriate low-energy theory for an individual SWNT 
is the Luttinger liquid (LL) model [H H HJ. The 
LL theory of SWNTs is usually formulated within the 
Abelian bosonization approach j^. With x = (x, t), 
where x is the spatial ID coordinate along the tube, and 
< r < 1/T is imaginary time, and corresponding in- 
tegration measure <ix = dxdr, the action for a single 
SWNT is [H[ll 

Sll = f dx V [(d T ip a /v a ) 2 + (d x ip a ) 2 } 

a—c±,s± J 

= J dx]T^ [(drOaM 2 + (d x 9 a ) 2 ] , (2.2) 
a 

which we take to be the same for every SWNT. Due to 
the electron spin and the additional K point degener- 
acy present in nanotubes 0], there are four channels, 
a — c+,c—,s+,s—, corresponding to the total/relative 
charge/spin modes [T^. l2*cj. with associated boson fields 
y a (x) and dual fields n (x) .21]. In the a = (c+,s— ) 
channels, the second (dual) formulation turns out to be 
more convenient, while the first line of Eq. 12. 2f) is more 
useful for a = (s+, c— ). The combined effect of Coulomb 
and phonon-mediated electron-electron interactions re- 
sults in the interaction parameter g c+ , where we assume 
g c+ > 1, reflecting effectively attractive interactions [Tl| . 
In the neutral channels, there are only very weak resid- 
ual interactions, and we therefore put g a ^c+ = 1 • Finally, 
the velocities v a in Eq. I|2.2|l are defined as v a — vp/g a , 
where vf = 8 x 10 5 m/sec is the Fermi velocity. 

Next we address the question which processes trigger 
the strongest superconducting fluctuations in a nanotube 
rope. This question has been addressed in Refs. 0, 0, 
Hi} . an d the conclusion of these studies is that Cooper 
pairs predominantly form on individual SWNTs rather 
than involving electrons on different SWNTs, see, e.g., 
the last section in Ref. ^3 f° r a detailed discussion. Fur- 
thermore, the dominant intra-tube fluctuations involve 
singlet (rather than triplet) Cooper pairs. The relevant 
order parameter for superconductivity is then given by 

E2 

0(x) = ^ 0'VV,<r, j 9(x)V>-r,-o-,- j 9(x), (2.3) 
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where ip r<r p denotes the electron field operator for a right - 
or left-moving electron (r — ±) with spin a = ± and K 
point degeneracy index (3 = ±. In bosonized language, 
this operator can be expressed as [2^] 



O = cos[v / 7T#c+] cos[v / 7rV3 c -] 

x cos[v / ttv 3 s+] cos[\/rr9 s -} — (cos «-» sin), 



(2.4) 



where we identify the UV cutoff necessary in the 
bosonization scheme with the graphite lattice constant, 
ao = 0.24 nm. In what follows, we use the short- 
hand notation <pj to label all four boson fields (p a (or 
their dual fields) corresponding to the jth SWNT, where 

J I v. 

The next step is to look at possible couplings among 
the individual SWNTs. In principle, three different pro- 
cesses should be taken into account, namely (i) direct 
Coulomb interactions, (ii) Josephson couplings, and (iii) 
single-electron hopping. The last process is strongly sup- 
pressed due to the generally different chirality of adja- 
cent tubes |T^ . and, in addition, for g c+ > 1, inter- 
S WNT Coulomb interactions are irrelevant ^| . Further- 
more, as discussed in detail in Ref. phonon-exchange 
mediated interactions between different SWNTs can al- 
ways be neglected against the intra-tube interactions. 
Therefore the most relevant mechanism is Josephson cou- 
pling between metallic SWNTs. These couplings define a 
Josephson matrix Aj k , which contains the amplitudes for 
Cooper pair hopping from the jth to the fcth SWNT. We 
put Ajj = 0, and hence A is a real, symmetric, and trace- 
less matrix. It therefore has only real eigenvalues A a , 
which we take in descending order, Ai > A2 > . . . > Ajv- 
Moreover, there is at least one positive and at least one 
negative eigenvalue. The largest eigenvalue Ai will be 
shown to determine the mean-field critical temperature 
Tj? below. The matrix A is then expressed in the corre- 
sponding orthonormal eigenbasis |a), 



(2.5) 



where (j\a) is the real orthogonal transformation from 
the basis of lattice points {\j}} to the basis {\a}} that 
diagonalizes A. Clearly, (j\a) = (ct\j). In what follows, 
we define ao such that A Q > for a < ao. 
The Euclidean action of the rope is then 



N 

3=1 jk 



(2.6) 



where Oj is the order parameter specified in Eq. I|2.4[) . 
The action l|2.6|l defines the model that is studied in the 
remainder of our paper. For studies of closely related 
models, see also Refs. [l3l 123) . 

In order to decouple the Josephson term in Eq. I|2.6|l . 
we employ a Hubbard-Stratonovich transformation. To 
that purpose, since the Josephson matrix has at least one 



negative eigenvalue, we first express A in its eigenbasis, 
see Eq. 1)2. 5p . The Josephson term in Eq. 12.2fl is then 
rewritten as 

J^O*A jk O k = £ sgn(A a )\A a \0* a O a , 

jk a 

where the order parameter in the |a) basis is 

O a =Y t (a\%)O i , 0* a = Y,0*(i\a). (2.7) 

i i 

By introducing a field A Q (x) for each Josephson eigen- 
mode |24j |. with (formally independent) complex con- 
jugate field A*, it is now possible to perform the 
Hubbard-Stratonovich transformation following the stan- 
dard procedure [25|. With integration measure VA = 
Y[ a 2?A* T>A a , the effective action entering the partition 
function Z = J VA exp(— S e g [A]) reads 

S off [A] = S [A] + J dx^A^A a , (2.8) 

where the action So [A] is formally defined via the re- 
maining path integral over the boson fields ipj, 



N 

So [A] = -In / n^e-£ 

J 3=1 



x e -;&E a c a (A;o t ,+o;Aj ) (29) 
with c a = 1 for a < ao, and c a — i otherwise. 



III. QUANTUM GINZBURG-LANDAU 
APPROACH 



A. Cumulant expansion 

Clearly, closed analytical evaluation of the path inte- 
gral in Eq. I|2.9|l is in general impossible. In order to 
make progress, approximations are necessary, and in the 
following we shall construct and analyze the Ginzburg- 
Landau (GL) action 0, |2^ for this problem. It turns 
out to be essential to take into account quantum fluc- 
tuations, i.e., the imaginary-time dependence of the or- 
der parameter field A q (x,t). In the standard (static) 
Ginzburg-Landau theory, such effects are ignored. 

The derivation of the GL action proceeds from a cumu- 
lant expansion of Eq. (|2.9|l up to quartic order in the A Q . 
This is a systematic expansion in the parameter | A|/27rT 
|25|. and by self-consistently computing this parameter, 
one can determine the regime of validity of GL theory. 
We stress that this expansion is not restricted to N ^> 1. 
In addition, for the long-wavelength low-energy regime 
of primary interest here, we are entitled to perform a 
gradient expansion. Using the single-chain correlation 
function G(xi 2 ) = (0(xi)0*(x 2 )) of the operator O in 
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Eq. (|2.4() with respect to the free boson action Sll, an d 
the connected four-point correlation function 

G( 4 )( Xl ,x 2 ,X3,x 4 ) - (0(x 1 )0(x 2 )0*(x 3 )0*(x 4 )) 

- (0(x 1 )O*(x 3 )){O(x 2 )O*(x 4 )) 

- (0(x 1 )0*(x 4 ))(0(x 2 )0*(x 3 )), 

the cumulant-plus-gradient expansion up to quartic order 
yields for the effective Lagrangian density 



with the real, symmetric, and positive definite matrix 



L\A] 



A a \ 2 + D\d T A a \ 



(3.1) 



a<ao 



(A-^IAJ 2 " 

B V f»i>«f A* A* A A 



where we use the notation 



i 

The temperature-dependent positive coefficients 
A,B,C,D are obtained as 



A 
B 
C 
D 



dxG(x), (3.2) 

_ 4 / dxiC?X2dx3 G c 4) ( x i: x 2,x 3 ,x 4 ), (3.3) 

~ /" dxx 2 G(x), (3.4) 

i | dxr 2 G(x). (3.5) 



Due to translation invariance, the integral for B does not 
depend on x 4 . Besides temperature, these coefficients 
basically depend only on the important LL interaction 
parameter g c+ . In particular, as it is discussed below, 
for g c+ > 1, the coefficient A grows as T is lowered. 
For static and uniform configurations, modes with a > 
«o never become critical. One can then safely integrate 
over these modes, which leads to a renormalization of 
the parameters governing the remaining modes. Such 
renormalization effects are however tiny, and thus are 
completely neglected in Eq. (|3.1|) . 

At this stage, it is useful to switch to an order param- 
eter field defined on the jth SWNT, 



(3.6) 



Q<ao 



After some algebra, the Lagrangian density (|3.1(l can be 
written as 

N 

L[A] =J2[c\d x A j \ 2 + D\d T Aj\ 2 + -B|Aj| 4 + 



3=1 



(A^-^IA/j +5>^ fc A fe) (3 . 7) 

jk 



v jk = E^WO^-VH* 



(3.8) 



a<Qo 



Notice that, strictly speaking, the fields Ai are not all in- 
dependent, because we have defined them from the sub- 
set of positive modes. The transformation in Eq. i|3.6|) 
is indeed not invertible. Nevertheless, in the following, 
we treat them as formally independent. This only affects 
the precise values of the Vij but does not qualitatively 
change our results. The expectation value of the order 
parameter field (|2.4I) can be expressed in terms of linear 
combinations of the fields Aj(x, t), and hence it is indeed 
justified to call Aj a proper "order parameter field". 

Equation l|3.7|l specifies the full GL action, taking into 
account quantum fluctuations and transverse modes for 
arbitrary number N of active SWNTs. In the limit 
N — > oo, and considering only static field configurations, 
results similar to those of Ref. [13j are recovered. In 
that limit the last term in Eq. I|3.7|l gives indeed the 
gradient term in the transverse direction, and one ob- 
tains the standard 3D GL Lagrangian. There is how- 
ever an important difference, namely the starting point 
of Ref. 13] is a model of Josephson-coupled ID supercon- 
ductors, whereas we start from an array of metallic chains 
with g c+ > 1, where the inter-chain Josephson coupling 
is crucial in stabilizing superconductivity. More similar 
to ours is the model investigated in Ref. |23j. However, 
in that paper, the metallic chains are assumed to have a 
spin gap, which is not the case for the SWNTs in a rope 
in the temperature range of interest. Furthermore, the 
main focus in Ref. p3| is the competition between charge 
density wave and superconducting instabilities, whereas 
in our case, as discussed above, the formation of a charge 
density wave is strongly suppressed by compositional dis- 
order, i.e., different chiralities of adjacent tubes, and we 
do not have to take the corresponding instability into 
account. 



B. Ginzburg-Landau coefficients 

In order to make quantitative predictions, it is neces- 
sary to compute the GL coefficients defined in Eas. H3.2|) - 
(|3.5() . While this is possible in principle for the full four- 
channel model 1)2. 2|l . here we will instead derive the coeffi- 
cients for a simpler model, where the K point degeneracy 
is neglected. This leads to an effective spin-1/2 Luttinger 
liquid action with interaction parameter g c (g s = 1) and 
velocity v c = vp / g c (v s = vp). Up to a prefactor of or- 
der unity, the respective results can be matched onto each 
other. This can be made explicit, e.g., for the coefficient 
A, where we get from the full action ill'. 21 

_ c_ ( ira T 

The proportionality constant c is found to differ from 
A/2ir 2 [see Eq. (|3.9[) below, which follows from the spin- 



1/2 description] only by a factor of order unity. In the 
simpler model neglecting the K point degeneracy, one 
then needs to take 



1 



9c 



which gives, for g c+ = 1.3, a value of g c w 1.1. This 
way, all exponents of the resulting power-law correlation 
functions (which are the physically relevant quantities) 
in the "reduced" model are the same as in the complete 
model, and only prefactors of order unity may be different 
for the respective GL coefficients. The bosonized order 
parameter 12.4fl in the simpler model is then given by 

O = cos[v / 27r cp s ] exp[i\/27r#c]- 

Using the finite-temperature correlation functions of the 
fields 8 C and ip s |2lj . 



(0 c (x)0 c (O)) = 



(tp s (-x.)<p s (0)) 



1 



2ng c 



In 



sinh 



ttT(x + iv c r) 



2n \TraoT 



sinh 



nT(x + ivft) 



v F 



and rescaling the integration variables x and r in 
Eqs. (j33-(E!3, explicit expressions follow in the form 



A(T) 
B(T) 
C{T) 
D{T) 



1 



2tt 2 vf 



(iraoT/v^ 1 - 1 !, 



(3.9) 



32n 4 v c v 2 F 



{ira T/v c ) g c ~ 3 C, 



■(naoT/v c ) 9 ° ^D. 



Dimcnsionless g c -dependent numbers A, B, C, D were de- 
fined as follows. With the notation z = (w,u) and 



dz = \ du I dw, 

o J -oo 



we have 



.4 
C 
D 



dz 



o 

/c(z)/ a (z)' 
„,2 



dz- 



where functions f c , s are introduced as 

f c (z) = |sinh(w; + m)| 1/9 % 
f s (z) = \smh(w/g c + iu)\. 



The coefficient of the quartic term in the GL functional 
is 



B = 



dzidz 2 dz 3 

/c(z 2 )/ c (zi 3 ) 

/ s (zi)/ s (z 23 ) 



4 / e (zi)/ c (z 23 ) 

/s(z2)/s(zi 3 ) / c (z 3 )/ c (zi 2 ) 

+ (1 4-» 2) + (1 <-> 3) 



/a (z 2 ) / a (z i3 ) / s (z 3 ) f s (z i 2 ) 



with Zjj = (u>i — Wj,m — Uj). The quantity 1? is eval- 
uated using the Monte Carlo method. For g c — 1, we 
first numerically reproduced the exact result B = 8ir 2 C 
with C = 77r£(3)/4 [25|. Numerical values can then be 
obtained for arbitrary g c . Numerical evaluation yields 
for g c 1.1 (corresponding to g c+ = 1.3) the following 
results: 



A ~ 17.4, B~ 392(1), C 



.15, 



D ~ 6.97. 

(3.10) 



C. Mean-field transition temperature 

Since in the rope only a modest number of transverse 
modes are present, a natural definition of the mean-field 
critical temperature Tj? is the temperature at which the 
mode corresponding to the largest eigenvalue of A be- 
comes critical. From Eq. (|3.7|l . this leads to the condi- 
tion A(T) = A^ 1 , and hence to the mean-field critical 
temperature 



T o = Vc_ I AAi 
c nao \ 2tt 2 v f 



9c/(3c-l) 



(3.11) 



which exhibits a dependence on the number N of active 
SWNTs in the rope through Ai. For large N, the eigen- 
value Ai saturates, and Eq. I|3.11|l approaches the bulk 
transition temperature. 

To provide concrete theoretical predictions for TJ? is 
difficult, since the Josephson matrix is in general un- 
known, and the results for T® very sensitively depend 
on Ai. Using estimates of Ref. [ljj and typical N as re- 
ported in Ref. Q, as an order-of- magnitude estimate, we 
find T° values around 0.1 to 1 K. When comparing to ex- 
perimental results, Ai can be inferred from the actually 
measured T c , which in turn provides values in reasonable 
agreement with theoretical expectations 0. 



D. 



Low-enerey theory: T < T c ° 



In what follows, we focus on temperatures T < T c °. 
Then it is useful to employ an amplitude-phase represen- 
tation of the order parameter field, 



A j( x ) = l A il( x )exp[^-(x) 



(3.12) 



where the amplitudes |A 3 -| are expected to be finite with 
a gap for fluctuations around their mean-field value. At 
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FIG. 1: Temperature dependence of A /2ttT versus T/T c ° 
for JV = 31 (open circles) and N = 253 (filled circles). 



^2^|A l ||A J |cos( 



with the Mooij-Schon velocity |26j . 



c s = v c \J C/D, 
and dimensionless phase stiffness parameters 
Hi =2^C|A J | 2 /c s . 



(3.14) 



(3.15) 



(3.16) 



At this stage, electromagnetic potentials can be coupled 
in via standard Peierls substitution rule > an d dissipa- 
tive effects due to the electromagnetic environment can 
be incorporated following Ref. (IT) ■ 



not too low temperatures, the GL action corresponding 
to Eq. I|3.7I) is accurate (see below), and the mean-field 
values follow from the saddle-point equations. Consider- 
ing only static and uniform field configurations, we find 
(j>i = cj>, where in principle also other (frustrated) con- 
figurations with <j)i — <f>j — ±7r could contribute. Such 
configurations presumably correspond to maxima of the 
free energy, and are ignored henceforth. The saddle-point 
equations then reduce to equations for the amplitudes 
alone, 



{A^ 1 -A)\A i \+2B\A i \ 3 = 0, (3.13) 



whose solution yields the transverse order parameter pro- 
file. Numerical study of Eq. (|3~T3|) using a standard 
Newton-Raphson root-finding algorithm then allows to 
extract the profile {|Aj|} for a given Josephson matrix 
Ay. We briefly discuss the solution of Eq. (|3.13|) for 
the idealized model of a rope as a trigonal lattice exclu- 
sively composed of N metallic SWNTs, where Ay = A 
for nearest neighbors and Ay = otherwise. For 

this model, Fig. ^ shows the resulting average ampli- 
tude Ao = Ylt \^i\/N as a function of temperature for 
X/vp = 0.1 and two values of N. Since Ao/2-7rT is the 
expansion parameter entering the construction of the GL 
functional, and it remains small down to T k /2, we 
conclude that the GL theory is self-consistently valid in 
a quantitative way down to such temperature scales. In 
our discussion below, GL theory turns out to be qualita- 
tively useful even down to T = 0. 



IV. ID ACTION AND QUANTUM PHASE 
SLIPS 



A. ID phase action 

Numerical evaluation of Eq. (|3.13(l shows that for T 
well below TJ?, transverse fluctuations are heavily sup- 
pressed. While this statement only applies to amplitude 
fluctuations, one can argue that also the transverse phase 
fluctuations are strongly suppressed. The basic argument 
relates to the scaling dimension [in the renormalization 
group (RG) sense] of the operator cos(0i — <j>j), which is 
essentially governed by the fij. For T well below T°, the 
become large, and the cosine operators get strongly 
relevant, locking the phases all together. In the low- 
temperature regime of main interest below, this argu- 
ment allows to substantially simplify Eq. (|3.14() . Then 
also no detailed knowledge about the Josephson matrix 
is required, because the only relevant information is es- 
sentially contained in T c °. 

Putting all phases <pj — <f>, we arrive at a standard 
(Gaussian) ID superconducting phase action |l4j . 



S 



JL 

2tt 



dxdr [^{d^f +C s {d x (j>f 



(4.1) 



with dimensionless rigidity fi = J^j^j^ see Eq. I|3.16|l . 
and c s as given in Eq. (|3.15|l . Assuming GL theory to 
work even down to T — for the moment, and neglecting 
the Vy-term in Eq. (|3.13(l . a simple analytical estimate 
follows in the form 



H(T) =Nv\l- (T/T^-V' 3 ' 



Fixing the amplitudes |Aj| at their mean-field val- where the number v is 



ues, and neglecting the massive amplitude fluctua- 
tions around these values, the Lagrangian follows from 



4ttA(CT>) 1/ 7B 



(4.2) 



(4.3) 
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The peculiar temperature dependence of the phase stiff- 
ness in Eq. 14. 2[) . reflecting the underlying LL physics 
of the individual SWNTs, is one of the main results 
of this paper. In the effective spin-f/2 description em- 
ployed here, using the numbers specified in Eq. (|3.10() 
for g c = 1.1 results in v « 4. Remarkably, at T = 0, 
Eq. (|4.2|) coincides, up to a prefactor of order unity, with 
the rigidity ft obtained from standard mean-field relations 



/i = TT 2 n s R 2 /2m*c s 



DN. 



With the density of condensed electrons n s and rope ra- 
dius R, this implies v ss vf/c s , which is of order unity. 
We therefore conclude that the GL prediction (|4.2[l for 
/J,{T) is robust and useful even outside its strict validity 
regime. 

The result l|4.2Jl for the stiffness is central for the fol- 
lowing discussion. The value we obtain for v, however, 
should not be taken as a very precise estimate. First, it 
can be affected by factors of order unity under a full four- 
channel calculation taking into account the K point de- 
generacy, as this affects each of the numbers in Eq. (|3.I()|1 
by a factor of order unity. Second, uncertainties in the 
parameter g c will also affect v by a factor of order unity. 
Moreover, based on the discussion in Ref. 0], one ex ~ 
pects on general grounds that intra-SWNT disorder and 
dissipative effects, both of which are not included in our 
model, will effectively lead to a decrease of the parame- 
ter v entering Eq. (|4.2[) . Therefore v is taken below as 
a fit parameter when comparing to experimental data. 
Since the number of active SWNTs N can be estimated 
from the residual resistance, and the transition temper- 
ature T c , see Eq. (|4.7|l below, can be determined from 
the experimentally observed transition temperature, v 
is basically the only free remaining parameter. Fits of 
our theoretical results to experimental data are then ex- 
pected to yield values for v around v w f . This is verified 
below in Sec. I VII 



B. Phase slips 

In the ID situation encountered here, superconductiv- 
ity can be destroyed by phase slips |14j. A phase slip 
(PS) can be visualized as a process in which fluctuations 
locally destroy the amplitude of the superconducting or- 
der parameter, which effectively disconnects the ID su- 
perconductor into two parts. Simultaneously, the phase, 
being defined only up to 2ir, is allowed to "slip" by 27r 
across the region where the amplitude vanishes. This 
process then leads to finite dissipation in the supercon- 
ducting wire via the Josephson effect. Depending on tem- 
perature, phase slips can be produced either by thermal 
or by quantum fluctuations. In the first case, which is 
commonly realized very near the critical temperature, we 
have a thermally activated phase slip (TAPS). At lower 
temperature, the quantum tunneling mechanism domi- 
nates, and one speaks of a quantum phase slip (QPS). 



For a textbook description of quantum phase slips, see 
Ref. (23. Below we demonstrate that in superconduct- 
ing ropes, only QPSs are expected to play a prominent 
role. 

A QPS is a topological vortex-like excitation of the su- 
perconducting phase field (/>(x, r) that solves the equa- 
tion of motion for the action (|4.I|I with a singularity 
at the core, where superconducting order is locally de- 
stroyed and a phase cannot be defined. Defining a ther- 
mal lengthscale as 



(4.4) 



for rope length L — > 00 and Lt — > 00, a QPS with core 
at (xi,Ti) and winding number ki = ±1 (higher winding 
numbers are irrelevant) is given by |2?| 



, t) = ki arctan 



c s (t - n) 



(4.5) 



where the finite L, Lt solution follows by conformal 
transformation [l^. The action of a QPS consists of 
two terms, one associated with the local loss of conden- 
sation energy, the core action S c , and the other with the 
vortex strain energy. While a detailed computation of S c 
requires a microscopic description of the dynamics inside 
the vortex core |l6j | , a simple qualitative argument is able 
to predict an order-of-magnitude estimate S c rs fi/2 |27| . 

This result allows us to assess the relative contribu- 
tion of the TAPS and QPS mechanisms. The produc- 
tion rate for the creation of one vortex is [Til 7qps ~ 

^ Cs exp(— 5 C ), where k is the core size. Within expo- 
nential accuracy, comparing this formula to the respec- 
tive standard TAPS rate expression [1J|, the crossover 
temperature from TAPS- to QPS-dominated behavior 
is Tp S = 2AF/Nv, with activation barrier AF. Us- 
ing results of Ref. 23] , we estimate the latter as AF = 
8V2R(g c )NT® /3, with dimensionless coefficient R(g c ) of 
order unity. Finally, this implies Tp S w Tj?. Since the 
true transition temperature T c < Xj?, see below, in the 
temperature regime T < T c , the influence of a TAPS can 
safely be neglected against the QPS. 

The generalization to many QPSs then leads to the 
standard picture of a Coulomb gas of charges ki = ±1, 
with fugacity y = e c , total charge zero, and log- 
arithmic interactions |25l |2?| . The partition function 
Z = ZqZv contains a regular factor Zq and the vor- 
tex contribution 



Z v = 



\ - y 



^ (n!) 2 

n=Q y ' 



(4.6) 

This model undergoes a Berezinski-Kosterlitz-Thouless 
transition driven by the nucleation of vortices, here cor- 
responding to a transition from a phase fj, > fi* , where 
QPSs are confined into neutral pairs and the rope forms a 
ID superconductor with finite phase stiffness and quasi- 
long-range order, to a phase fi < fj,* where QPSs prolif- 
erate. In that phase, vortices are deconfined and destroy 



8 




FIG. 2: Temperature-dependent resistance R(T < T c ) pre- 
dicted by Eq. 15.41 for v = 1 and different N. The smaller 
is N, the broader is the transition. From the leftmost to the 
rightmost curve, N = 4, 7, 19, 37, 61, 91, 127, 169, 217. 



the phase stiffness, thereby producing normal behavior, 
where "normal" does of course not imply Fermi- liquid be- 
havior. The phase boundary is located at fi* = 2 + Airy ~ 
2. The true transition temperature T c is therefore not the 
mean- field transition temperature T c °, but follows from 
the condition fi(T c ) = fx*. Putting fi* = 2, Eq. (|4"2|) 
yields 



T c /T c ° = [1 - 2/Nu 



ga/iga-l) 



(4.7) 



This T c depression is quite sizeable for N < 100. To give 
concrete numbers, taking v = 1, for N = 25, 50, and 100, 
the ratio T c /T c ° equals 0.40, 0.63, and 0.80, respectively. 
QPSs also have an important and observable effect in the 
superconducting regime, as will be discussed in the next 
section. 



V. RESISTANCE BELOW T c 

A phase slip produces finite dissipation through the 
Josephson effect, and therefore introduces a finite resis- 
tance even in the superconducting state, T < T c . The 
QPS-induced linear resistance R(T) = V/I for T < T c 
can be computed perturbatively in the QPS fugacity y 
[l5j . For that purpose, we imagine that one imposes a 
small current / to flow through the rope. The presence 
of QPSs implies that a voltage drop V occurs, which is 
related to the average change in phase, 



V 



2e 



-[r(/)-r(-/)], 



where T(±I) is the rate for a phase slip by ±27r |15| . 
This rate can be obtained following Langer [2^ as the 
imaginary part acquired by the free energy F(I) under 
an appropriate analytic continuation, 



F(±J) = -21mF(±I). 



(5.1) 



We only consider the contribution of a single pair of 
QPSs, i.e., compute R(T) to second order in y. Expand- 
ing Eq. H4.6(l to order y 2 , the free energy at this order 
reads 



F 



Ly 2 c 2 



l/T j-L/2 

dr / dxe eT ~ 2t " 9E( - x ' T \ (5.2) 

i' J-L/2 



where the vortex- vortex interaction gE(x,r) only de- 
pends on relative coordinates, and e = irhl/e. The con- 
tribution Fq to the free energy due to regular configu- 
rations can be dropped, because it does not acquire an 
imaginary part under the analytic continuation. We now 
perform the analytic continuation r — > it, resulting in 
[see Ref. [H for details] 



ImF = 



T 2 2 

2k 4 



L/2 



dx 



L/2 



where g(x, t) = gE^x, r — > it). The rate T(e) then follows 
for L, Lt S> k but arbitrary L/ Lt in the form 



r(e) 



where 



c 2 s Ly 2 



L/2 



dx 



L/2 



g(t) = In [(L t /k) sinh(7rT|t|)] + i(7r/2)sgn(t). 

Analyticity of gE(x,r) in the strip < r < l/T also 
leads to the standard detailed balance relation [23 , 

r(-e) = e-^ T T(e). 

In order to explicitly evaluate the rate T(e) for arbitrary 
L/Lt, we now replace the boundaries for the x-integral 
by a soft exponential cutoff, switch to integration vari- 
ables t' = t — x/c s and t" = t+x/c s , and use the auxiliary 
relation 



-c s \t"-t'\/L _ _£f_ 

nL 



ds- 



-is(t" -t') 

+ {c s /L) 2 



The t',t" time integrations then decouple, and it is 
straightforward to carry them out. Finally some alge- 
bra yields the linear resistance in the form 



R(T) = / iryTQi/2) 

dll : 



2k V K J 



r(/i/2 + iuL T /2L) 



Tin/2) 



in units of the resistance quantum Rq defined in 
Eq. f2~TJ) . 

For L/Lt 3> 1, the it-integral approaches unity, and 
hence R cx T 2m ~ 3 , while for L/Lt <§C 1, dimensional scal- 
ing arguments give R cx T 2 ^- 2 . In Refs. 0,0, typical 
lengths were Lwl /im, which puts one into the crossover 
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FIG. 3: Temperature dependence of the resistance below 
T c for superconducting rope R2 experimentally studied in 
Ref. Q. Open squares denote experimental data (with sub- 
tracted residual resistance), the curve is the theoretical result. 

regime Lt ~ L. While the quoted power laws have al- 
ready been reported for diffusive wires 0] , Eq. (|5.4|) de- 
scribes the full crossover for arbitrary L/Lt, and applies 
to strongly correlated ladder compounds such as nan- 
otube ropes. It predicts that the transition gets signifi- 
cantly broader upon decreasing the number of tubes in 
the rope. This is shown in Fig. |3 where the theoreti- 
cal results for the resistance is plotted for various N at 
v = 1. Note that Eq. I|5.4|l is a perturbative result in the 
fugacity, and it is expected to break down close to T c , see 
below. In the next section we directly compare Eq. (15. 4|) 
to experimental data obtained by Kasumov et al. jj]. 

VI. COMPARISON TO EXPERIMENTAL DATA 

Here we discuss how the prediction for the 
temperature-dependent resistance R(T) below T c as 
given in Eq. I|5.4|) compares to the experimental results 
for R(T) published in Ref. More aspects of this com- 
parison will be given elsewhere 30]. The experimental 
data in Ref. Q were obtained from two-terminal mea- 
surements of ropes suspended between normal electrodes. 
Due to the presence of the contacts, the residual resis- 
tance survives down to T = even when the rope 
exhibits a superconducting transition. Extrapolation of 
experimental results for R(T) yields R Cl which then al- 
lows to infer the number N in the respective sample 
from Eq. (|2.1|) . This resistance R c has to be subtracted 
from experimental data to allow for a comparison with 
Eq. I|5.4I) . where no contact resistance is taken into ac- 
count. 

In Figs. |3| and experimental resistance curves (af- 
ter this subtraction) for the samples named R2 and R4 
in Ref. Q are plotted versus the prediction of Eq. I|5.4II . 
For sample R2, we find R c = 74 ft corresponding to 
Nr2 = 87, while for sample i?4, the subtracted resis- 



tance is R c = 150 ft, leading to Nr4 = 43. We then take 
these N values when computing the respective theoreti- 
cal curves. The experimentally determined temperature 
T* locates the onset of the transition 4], and is identified 
with the true transition temperature T c in Eq. 14. 7f) . It 
is therefore also not a free parameter. Note that thereby 
the eigenvalue Ai of the Josephson matrix has been de- 
termined. In the absence of detailed knowledge about 
the Josephson matrix, it is fortunate that our result for 
R{T)/R{T C ) following from Eqs. JOJ| and 63 docs not 
require more information about A besides the largest 
eigenvalue. Given the estimate g c + = 1.3 11], the com- 
parison of Eq. I|5.4I) to experimental data then allows only 
one free fit parameter, namely v. According to our dis- 
cussion in Sec. IIV Al the fit is expected to yield values 

The best fit to the low-temperature experimental 
curves for R(T) yields v = 0.75 for sample R2, see Fig. |3I 
and v = 0.16 for sample i?4, respectively. The agreement 
between experiment and theory is excellent for sample 
R2. For sample i?4, the optimal v is slightly smaller 
than expected, which indicates that dissipative processes 
may be more important in that sample. Nevertheless, 
for both samples, the low-temperature resistance agrees 
quite well, with only one free fit parameter that is found 
to be of order unity as expected. Whereas the agree- 
ment between theoretical and experimental curves ap- 
pears then quite satisfactory in the low-temperature re- 
gion, our predictions clearly deviate in the region near 
T c . This is not surprising, because our expression for 
R(T) in Eq. I|5.4|l is perturbative in the QPS fugacity. It 
is then expected to break down close to T c , where QPSs 
proliferate and the approximation of a very dilute gas of 
QPS pairs, on which our calculation is based, is not valid 
anymore. As a consequence, also the saturation observed 
experimentally above T* is not captured. 

We note that it is an interesting challenge to compute 
the finite resistance in the normal phase at T c < T < Tj?, 
where the saturation should be caused by QPS and TAPS 
proliferation. For temperatures T > T®, superconduct- 
ing correlations can be neglected, and the resistance 
should then be dominated by phonon backscattering and 
disorder effects. Nevertheless, we believe that the agree- 
ment between the theoretical resistance result (|5.4f) and 
experimental data at low temperatures shown in Figs. |3| 
and^J given the complexity of this system, is rather sat- 
isfactory. More importantly, this comparison provides 
strong evidence for the presence of quantum phase slips 
in superconducting nanotube ropes. 



VII. CONCLUSIONS 

According to our discussion above, the intrinsic super- 
conductivity observed in ropes of SNWTs 0, Q repre- 
sents a remarkable phenomenon, where it has been pos- 
sible to experimentally probe the extreme ID limit of a 
superconductor. In this paper, we have formulated a the- 
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FIG. 4: Same as Fig. [3] but for sample R4 experimentally 
studied in Ref. Q. 



Based on the resulting low-energy action for the phase 
fluctuations, we have shown that quantum phase slips 
produce a depression of the critical temperature. More 
importantly, the temperature dependence of the linear 
resistance experimentally observed below the transition 
temperature can be accounted for by considering the un- 
derlying LL physics and the effect of quantum phase 
slips. Despite some admittedly crude approximations, 
like the neglect of intra-tube disorder and dissipation ef- 
fects inside the vortex core, the comparison of experimen- 
tal curves and theoretical predictions, in particular in the 
low-temperature region, strongly suggests that the resis- 
tive process is indeed dominated by quantum phase slips. 
Our theory also suggests that, if repulsive Coulomb inter- 
actions can be efficiently screened off, superconductivity 
may survive down to only very few transverse channels 
in clean nanotube ropes. 



ory for this phenomenon, based on a model of metallic 
SWNTs with attractive intra-tube interactions and ar- 
bitrary inter-tube Josephson couplings. The analysis of 
this model leads to an effective Ginzburg-Landau action, 
whose coefficients can be expressed in terms of parame- 
ters entering the microscopic description of the rope. In 
order to get the correct low-energy dynamics, it is crucial 
to include quantum fluctuations of the order parameter. 
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